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ABSTRACT 

We describe a correlation between the mass M, of a galaxy's central black hole and the luminosity-weighted 
line-of-sight velocity dispersion a e within the half-light radius. The result is based on a sample of 26 galaxies, in- 
cluding 13 galaxies with new determinations of black hole masses from Hubble Space Telescope measurements of 
stellar kinematics. The best-fit correlation is M. = 1.2(±0.2) x lO s M Q (a e /200km s -!) 3 - 75 (± - 3 ) over almost three 
orders of magnitude in M.; the scatter in M, at fixed <j e is only 0.30 dex and most of this is due to observational 
errors. The M,-a e relation is of interest not only for its strong predictive power but also because it implies that 
central black hole mass is constrained by and closely related to properties of the host galaxy's bulge. 

Subject headings: galaxies: nuclei — galaxies: statistics — galaxies: general 



1. INTRODUCTION 

Massive black holes at the centers of galaxies are now rec- 
ognized as a normal, perhaps ubiquitous, component of ellip- 
tical galaxies and spiral galaxy bulges. The early evidence 
is summarized in Kormendy & Richstone (1995) and refer- 
ences therein. Using a heterogeneous set of 36 galaxies, mostly 
with ground-based spectroscopy and space-based photome- 
try, Magorrian et al. (1998) argue that all hot galaxy compo- 
nents (ellipticals and spiral galaxy bulges) contain central black 
holes. Kormendy (1993), Kormendy & Richstone (1995), and 
Magorrian et al. (1998) also find that black hole mass M m is 
proportional to galaxy mass or luminosity, although with large 
scatter (see left panel of Fig. 2). Richstone et al. (1998) out- 
line a plausible physical framework to discuss the connections 
between current black hole mass, galaxy formation, and quasar 
evolution. The set of galaxies with reliable black hole masses 
is growing rapidly, mostly through an aggressive campaign of 
Hubble Space Telescope (HST) observations, and this set is now 
large enough to investigate the dependence of M. on galaxy 
properties, and thus to provide clues to the role of central black 
holes in galaxy formation and evolution. In this Letter, we 
present a new correlation between line-of-sight velocity disper- 
sion and black hole mass that has very little intrinsic scatter 
(probably less than 40% in M.). 

2. THE SAMPLE 

Our sample is restricted to galaxies whose black hole masses 
and line-of-sight velocity dispersions are well determined. We 
use only two galaxies with maser masses (NGC4258 and 
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NGC 1068), since only these have well measured dispersions. 
We also use six galaxies with black hole masses determined 
from gas kinematics. Among the galaxies with black hole 
masses determined from stellar kinematics, we limit our sample 
to those that have three-integral models with HST spectroscopy 
(16 galaxies), plus our Galaxy (Ghez et al. 1998; Genzel et al. 
2000) and M31 (Dressier & Richstone 1988, Kormendy 1988). 
We do not use 21 galaxies from the sample of Magorrian et al. 
(1998) that have only ground-based kinematic data and two- 
integral dynamical models. Galaxies with black hole detec- 
tions based on isotropic or spherical models (e.g., NGC 31 15: 
Kormendy et al. 1996a, NGC 4594: Kormendy et al. 1996b, 
NGC4486B: Kormendy et al. 1997) are also excluded until 
they are re-analyzed with three-integral models. 

The majority of black hole mass estimates in the present sam- 
ple (14 out of 26) come from galaxies with HST spectra and 
preliminary three-integral axisymmetric dynamical models de- 
scribed by Gebhardt et al. (2000a, 2000b) and Richstone et al. 
(2000). Since all of these galaxies have been observed, ana- 
lyzed, and modeled with the same procedures, we expect that 
the scatter due to systematic errors will be smaller than in a 
more heterogeneous sample. 

Table 1 presents the data as follows: galaxy name (Col- 
umn 1), type (Col. 2), black hole mass and 68% uncertainty 
(Col. 3), integrated line-of-sight velocity dispersion as defined 
below (Col. 4), distance in Mpc (Col. 5), and source (Col. 6). 
The velocity dispersions are from heterogeneous sources with 
uncertain errors, but probably most are accurate to within 
±5%. Most distances come from surface brightness fluctua- 
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tions (Tonry et al. 2000); otherwise we use Virgo infall models 
to infer the distance (Dekel 2000, private communication), all 
scaled to a Hubble constant of 80 km s -1 Mpc -1 . 



Table 1 
Galaxy Sample 



Galaxy 


Type 


M BH (Low.High) 




Dist 


Sources 








Stellar Dynamical Estimates 






Galaxy 


Sbc 


p 


2.6xl0 6 


(2.4,2.8) 


75 


0.008 


GeOO 


M31 


Sb 


P 


3.5xl0 7 


(1.0,6.0) 


160 


0.8 


K88,v94 


M32 


E2 


P 


3.7xl0 6 


(2.4,4.4) 


75 


0.8 


v98,v94 


N821 


E4 


P 


5.0xl0 7 


(3.0,7.0) 


209 


24.1 


GOOb 


N1023 


SO 


P 


3.9xl0 7 


(2.8,4.8) 


205 


11.4 


BOO 


N2778 


E2 


P 


2.0 xlO 7 


(0.7,3.6) 


175 


22.9 


GOOb 


N3377 


E5 


P 


1.0x10 s 


(0.6,2.5) 


145 


11.2 


G00b,K98 


N3379 


El 


c 


1. Ox 10 s 


(0.5,1.6) 


206 


10.4 


GOO a 


N3384 


SO 


P 


1.8 xlO 7 


(0.9,2.5) 


143 


1 1.6 


GOOb 


IN jOUo 


E2 


c 


1.1x10 s 


(0.8,2.5) 


182 


ZJ.U 


OUU D 


N4291 


E2 


c 


1.5x10 s 


(0.8,4.5) 


242 


26.2 


GOOb 


N4342 


SO 


p 


3.0x10 s 


(2.0,4.7) 


225 


15.3 


C99 


N4473 


E5 


c 


1.0 x10 s 


(0.4,1.8) 


190 


15.7 


GOOb 


N4564 


E3 


p 


5.7 xlO 7 


(4.0,7.0) 


162 


15.0 


GOOb 


N4649 


El 


c 


2.0x10" 


(1.0,2.5) 


375 


16.8 


GOOb 


N4697 


E4 


p 


1.2x10 s 


(0.8,1.3) 


177 


11.7 


GOOb 


N5845 


E 


p 


3.2x10 s 


(2.5,5.0) 


234 


25.9 


GOOb 


N7457 


SO 


p 


3.4xl0 6 


(1.7,6.0) 


67 


13.2 


GOOb 








Gas Dynamical Estimates 






M87 


E0 


c 


2.5 xlO 9 


(1.8,3.2) 


375 


16.1 


H94,v94b 


N4261 


E2 


c 


5.4x10 s 


(4.3,6.3) 


315 


31.6 


F96,v90 


N4374 


El 


c 


1.8xl0 9 


(0.9,2.8) 


296 


18.4 


B98 


N6251 


E2 


c 


6.0 x10 s 


(2.0,8.0) 


290 


106.0 


F99.S90 


N7052 


E4 


p 


3.3x10 s 


(2.0,5.6) 


266 


58.7 


v98b 


IC1459 


E3 


c 


3.5 x10 s 


(1.4,4.8) 


323 


29.2 


K00 








Maser Dynamical Estimates 






N1068 


Sb 


p 


1.7xl0 7 


(1.0,3.0) 


151 


15.0 


G96,N95 


N4258 


Sbc 


p 


4.2xl0 7 


(4.0,4.4) 


120 


7.2 


M95.H98 



NOTE. — The type refers to Hubble type, and either power-law (p) or 
core (c) galaxy (see Faber et al. 1997). The black hole masses are in solar 
masses, dispersions are in km s~ ] , and distances are in Mpc. The first source 
listing is for the black hole mass, and the second is for the dispersion; if only 
one is given then both values are from the same source. The source citations 
are in the references. 

3. ESTIMATING THE VELOCITY DISPERSION 

We devote particular care to finding a suitable line-of-sight 
dispersion. The traditional estimates of galaxy velocity disper- 
sions come from the central regions (Faber et al. 1989; J0r- 
gensen & Franx 1994), but these are affected by M. in some 
galaxies. We focus on the aperture dispersion (the luminosity- 
weighted line-of-sight dispersion inside a radius R), since it has 
a higher signal-to-noise and is less sensitive to the details of the 
distribution of orbits (Richstone & Tremaine 1984) than is the 
central dispersion. 

For the 13 galaxies studied by Gebhardt et al. (2000b) and 
Bower et al. (2000), we have extensive ground-based obser- 
vations of the line-of-sight dispersion as a function of radius. 
Figure 1 presents the aperture dispersion for these galaxies as 
a function of radius, normalized to the effective or half-light 
radius R e . The largest variations in the normalized dispersion 
profile occur in the center (R < Q.5R e ) and the outer regions 
(R > 2R e ). The variations among the dispersion profiles in 
the region 0.5R e < R < 2R e are small, generally less than 5%. 
This result follows because galaxies do not show dramatic ra- 
dial variations in dispersion, especially when integrated over 
an aperture; the same is generally true of theoretical models of 
spherical galaxies. In what follows, we measure "dispersion" 
by the line-of-sight aperture dispersion a e within R e . 




R/Re 

FIG. 1. — Luminosity-weighted line-of-sight dispersion within an aperture 
of radius R, normalized to its value at the effective radius R e . Each of the 13 
lines represents a galaxy from Gebhardt et al. (2000b) and Bower et al. (2000). 

The aperture dispersion includes a contribution from the ro- 
tation (i.e., we measure the rms velocity relative to the systemic 
velocity, not relative to the local mean velocity). The measured 
dispersion depends on the inclination of the galaxy. The values 
that we quote should be regarded as dispersions for galaxies 
as viewed edge-on, since several of the galaxies in the sample 
with the strongest rotation are nearly edge-on (the Galaxy, M3 1 , 
NGC 1023, NGC3377, NGC3384, and NGC4342). For the 
rest, the inclination is not well determined; however, the change 
in dispersion with viewing angle for ellipticals as judged from 
models is always < 20% and on average is less than 5% (as 
calculated from van der Marel 1991). 

4. RESULTS 

Figure 2 plots the black hole mass M. versus bulge luminos- 
ity and the aperture dispersion a e at R e . The bulge luminosity 
correlation will be discussed in Kormendy et al. (2000b). The 
correlation with dipersion is extremely strong, with a correla- 
tion coefficient of 0.93. By a variety of tests, it is much more 
significant than the 99% level. Furthermore, the correlation re- 
mains present in various subsamples — e.g., only galaxies from 
Gebhardt et al. (2000b), or only galaxies with gas dynamical or 
maser mass estimates. The correlation is robust. 

The line in the right panel Figure 2 is a fit to the data assum- 
ing that errors in dispersion measurements are zero and that 
errors in logM. are the same for each galaxy. The best- fit line 
is 

M - = L2Xl ° 8MQ (2^k^) 3 ' 75 - (1) 



We have also measured the relation using the seven galaxies 
with the best measured masses and dispersions (our Galaxy, 
M32, N4261, N4486, N4564, N4697, and N7052), and find the 
same fit within the errors. The 68% confidence limits for the 
best-fit linear relation come from standard least-squares fitting. 
The uncertainty, based on Monte Carlo simulations, is ±0.3 in 
the exponent and 0.057 dex in the normalization at 200 km s" 1 . 



BH Mass vs. Velocity Dispersion 
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FIG. 2. — Black hole mass versus bulge luminosity (left panel) and the luminosity-weighted aperture dispersion within the effective radius (right panel). There 
are 26 points in the dispersion plot; 13 are new detections from stellar kinematics (Gebhardt et al. 2000b, Bower et at 2000). Green squares denote galaxies with 
maser detections, red triangles come from gas kinematics, and blue circles are from stellar kinematics. Solid and dotted lines are the best-fit correlations and their 
68% confidence bands. 



The measured scatter is 0.30 dex in black hole mass at 
fixed dispersion. For the scatter estimate we use the biweight 
(Beers et al. 1990), and for this sample the biweight value is 
close to that of the standard deviation. Given the likely mea- 
surement errors, the intrinsic scatter is probably < 0.15 dex. In 
fact, the scatter from the seven galaxies with the best-measured 
quantities is 0.14 dex in black hole mass. More accurate black 
hole masses are required to better constrain the intrinsic scatter. 
With present data the intrinsic scatter is consistent with zero. 

It is possible that the tight correlation between black hole 
mass and dispersion arises in part because the sample that we 
analyze consists mainly of galaxies in HST observing programs 
that have been chosen to look "normal" (e.g., no unusual mor- 
phology, active nuclei, disturbed gas kinematics). It will be 
important to test the correlation using a broader sample. 

The galaxies that have the most leverage on the slope of the 
best-fit correlation are those at the extremes, both low and high 
mass. We discuss each of these, but we note that, when the 
three smallest and two largest black holes are excluded, the 
best-fit slope still lies within the original error bars. The three 
galaxies with the smallest black hole masses are the Milky Way, 
M32, and NGC 7457. The masses for our Galaxy (Ghez et al. 
1998; Genzel et al. 2000) and M32 (van der Marel et al. 1998) 
are among the most accurate of any in the sample. The mass 
for NGC 7457 (discussed in Gebhardt et al. 2000b) is diffi- 
cult to measure because there is a point source in the nucleus 
(Lauer et al. 1991); it may be overestimated. Measuring the 
line-of-sight dispersions for NGC 7457 and M32 is straightfor- 
ward (Pinkney et al. 2000b; van der Marel et al. 1998) since in 
both galaxies the dispersions are nearly independent of radius. 
The situation for our Galaxy is more difficult since we have 
to use the near-infrared surface-brightness profile (Kent et al. 
1991) and individual stellar velocities (Kent 1991; Genzel et al. 
2000) to determine the aperture dispersion. The effective radius 
of the bulge is uncertain due to disk contamination, but choos- 
ing values between 50 and 500 arcseconds provides dispersions 
between 85 and 65 km s" 1 . We adopt 75 km s" 1 . 

The two galaxies with the largest black holes are M87 



(NGC 4486) and NGC 4649. The black hole mass in M87 ap- 
pears well determined. The mass of the black hole in NGC 4649 
has substantial uncertainties since its nuclear velocity disper- 
sion is so large that absorption lines in the STIS spectrum are 
washed out (Gebhardt et al. 2000b). The effective velocity dis- 
persions in both of these galaxies have been measured by sev- 
eral groups with good agreement. 

We have looked for ways to reduce the scatter in the relation 
still further, but without much success. Including the effective 
radius as an additional parameter does provide a small improve- 
ment but only if the dependence is very weak, M. cx ^~ 01 at 
fixed a e . We also examined the correlation of M, with aperture 
dispersions within various radii, and find a similar correlation 
but with slightly larger scatter — by 0.03 dex. 

5. DISCUSSION 

A preliminary version of these results was presented at 
the June 2000 American Astronomical Society meeting in 
Rochester (Kormendy et al. 2000), where we learned of a sim- 
ilar analysis by Ferrarese & Merritt (2000). Their result, based 
on 12 objects, has a steeper slope (M. oc a 5 2 ). This differ- 
ence may be due to different dispersions, black hole masses, or 
distances for the same objects. Our sample size is larger (26 
compared to 12 galaxies), and many of our black hole masses 
come from as-yet-unpublished high-quality HST data and anal- 
yses. For their nearest galaxies (the Galaxy and M32) it is likely 
that the central dispersions used by Ferrarese & Merritt (2000) 
are enhanced by the black hole, which in turn creates a steeper 
slope. If we restrict our sample to the 1 1 galaxies common to 
their sample and ours but use our black hole masses and disper- 
sions, we find the same slope as presented here (M, oc crj 75 ). 
Thus, we are confident that our result is robust. 

The implications of our results can be discussed in the con- 
text of the fundamental plane for galaxies. Roughly speak- 
ing, elliptical galaxies and bulges can be described by three 
parameters: the effective radius R e , the total luminosity L, 
and the luminosity-weighted velocity dispersion a. If these 
galaxies have similar luminosity and mass distributions and if 



4 



Gebhardt et al. 



the mass-to-light ratio is a well behaved function of R e and 
L, then the virial theorem implies that they occupy a two- 
dimensional manifold in the three-dimensional space with coor- 
dinates (logL,loger,log/? c ); this is the so-called "fundamental 
plane" (Dressier et al. 1987; Djorgovski & Davis 1987; Faber 
et al. 1987). Any two of these variables can be used to predict 
the value of the third. 

Now let us examine the location of elliptical galaxies 
and bulges in the four-dimensional space with coordinates 
(logM. ,logL,loger,log/? e ). The correlation of M. with a that 
we describe here implies that (i) galaxies are still restricted 
to a two-dimensional manifold — a fundamental plane — in this 
four-dimensional space, and (ii) when projected onto the 
(log M m , log a) plane, the fundamental plane is viewed nearly 
edge-on — in other words, contours of constant M. on the fun- 
damental plane are parallel to contours of constant a. It is re- 
markable that the scatter normal to the plane, 0.05 dex in three 
dimensions, is increased to only 0.08 dex in four dimensions. 

Most earlier discussions of the demography of central black 
holes focused on the correlation of black hole mass with galaxy 
luminosity, which can be regarded as a projection of the four- 
dimensional fundamental plane onto the (logM.,logL) plane. 
We now understand that most of the substantial scatter in that 
relation — 0.6 dex in logM. — is not a reflection of stochastic 
processes that controlled the growth of black holes but arises 
instead simply because in this projection the fundamental plane 
is not viewed edge-on. Confirmation of this point is provided 
by the observation that those galaxies which are outliers in the 
(M. , L) relation are not outliers in the (M. , a e ) relation. 

The tight correlation between black hole mass and velocity 
dispersion strongly suggests a causal connection between the 
formation and evolution of the black hole and the bulge. How- 
ever, the nature of this connection remains obscure (Haehnelt & 
Rees 1993; Haiman & Loeb 1998; Silk & Rees 1998; Kauffman 
& Haehnelt 2000). It is natural to assume that bulges, black 



holes, and quasars formed, grew, or turned on as parts of the 
same process, in part because the collapse or merger of bulges 
might provide a rich fuel supply to a centrally located black 
hole. The finding that black hole mass correlates more closely 
with dispersion than total luminosity or other global properties 
probably reflects the fact that the aperture dispersion is less sen- 
sitive to the properties of the outer galaxy. 

We have shown that black hole mass is tightly coupled to ve- 
locity dispersion of the host galaxy over three orders of mag- 
nitude. It remains to be seen whether this correlation ap- 
plies at still larger or smaller black hole masses. It should 
be straightforward to explore smaller masses, but examining 
higher-dispersion galaxies (e.g., cD galaxies) will be challeng- 
ing because of their low surface brightnesses and large disper- 
sions. We may even speculate whether the correlation extends 
to stellar systems with dispersions as low as dwarf spheroidal 
galaxies and globular clusters, although in other respects these 
objects are quite different from bulges and elliptical galaxies 
(Kormendy 1985). A typical dwarf spheroidal galaxy or mas- 
sive globular cluster has a dispersion of about 10 km s" 1 ; if 
the correlation applies, it suggests a black hole mass of about 
2 x 10 3 M Q . It is intriguing that Gebhardt et al. (2000c) mea- 
sure an increase in the mass-to-light ratio of Ml 5 near its center, 
consistent with a central mass of 2500 M . 
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This work was supported by HST grants to the Nukers, GO- 
02600.01-87 A, G06099, and G07388, and by NASA grant G5- 
8232. A.V.F. acknowledges NASA grant NAG5-3556. K.G. 
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